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TECHNICAL NOTE NO. I559 

BUCKLING STRESSES OF CLAMPED RECTANGULAR 
FLAT PLATES IN SHEAR 

By Bernard Budiansky and Robert W. Connor 
SUMMARY 


By consideration of anti symmetrical, as well as symmetrical, 
■buckling configurations, the theoretical shear "buckling stresses of 
clamped rectangular flat plates are evaluated more correctly than in 
previous work. The results given, which represent the average of upper— 
and lower— limit solutions obtained by the Lagrapgian multiplier method, 

are within 1^ percent of the true buckling stresses. 


INTRODUCTION 


The theoretical buckling stress In shear of a clamped rectangular 
plate has been found exactly only for the case of an infinitely long 
plate (reference 1); for plates of finite aspect ratio, approximate 
solutions of uncertain accuracy have been presented by several investi- 
gators. In reference 2, Cox gives buckling— stress coefficients computed 
by the Eayleigh-Ritz energy method for several, aspect ratios, but no 
details of solution are Included. Iguchi (reference 3 V obtains approxi- 
mate results by means of a series method that, as Smith points out in 
reference 4, does not provide definite Information as to whether the 
results obtained are too high or too low. Smith then obtains upper 
limits to the true buckling stress by using Iguchi* s deflection function 
in the Rayleigh-Ritz method. The numerical results published by these 
investigators differ by as much as 10 percent at some values of , plate 
aspect ratio. Furthermore, these results were all apparently based on 
the assumption that, for each aspect ratio, the critical stress 
corresponds to a buckling pattern symmetrical (rather than anti symmetrical) 
about the plate midpoint. 

The present paper gives theoretical buckling— stress coefficients 
computed from an analysis (reference 5 ) by the Lagrangian multiplier 
method. Both upper and lower limits to the true buckling stress were 
calculated, so that the maximum error in the final results is definitely 
known. Furthermore, in the present paper the assumption is not made 
that the symmetrical buckling pattern governs at all aspect ratios. 
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In reference 6 It is shown that, in shear "buckling of a rectangular plate 
simply supported at the edges, an anti symmetrical "buckling pattern Is 
critical at certain aspect ratios. The ant I symmetrical pattern was 
accordingly included in the present-investigation, and the results 
indicate a range of aspect ratios in which the lowest buckling stress 
does occur with a configuration of antisymmetry. 


SYMBOIS 


a length of plate 

b width of plate 

3 plate aspect ratio (a/b) 

t thickness of plate 

TE Young* s modulus for material 

U Poisson's ratio for material 

D flexural stiffness of plate 

t critical shear stress 


Et 3 


. 12(1 - 


critical shear stress coefficient in the formula t = k. 


n ?D 
8 b2 t 


RESULTS AND DISCUSSION 


The 
edges is 


critical shear stress 
given by the equation 


for a rectangular flat plate with clamped 


T 


= k c 



The solid curve in figure 1 gives the values of the shear stress 
coefficient k B for values of aspect ratio 3 from 1 to 3 . As shown 

in figure 1 , the solid curve consists of parts of two distinct curves, 
one of which corresponds to symmetrical, buckling, the other to anti- 
symmetrical buckling, about the plate midpoint^— The governing 
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configuration at any aspect ratio, given "by the solid line, is that 
vhich yields the lover "buckling stress. 

The values of k B plotted in figure 1 represent the average of 
upper— and lover-limit solutions obtained from the analysis of 
reference 5 . (See table 1.) In most cases, the lover-limit results vere 
obtained from eleventh-order stability determinants, and the upper— limit 
results from ninthr-order determinants. (See reference 5*) The evaluation 
of many terms in the determinants vas considerably simplified by means 
of the confutation aids discussed in the appendix to the present paper. 

As is seen from the data of table 1, the final results obtained must be 

within 1 ^ percent of the true buckling stress coefficients. 

It is to be expected that for values of 0 between 3 and » the 
buckling stresses corresponding to symmetrical and anti symmetrical 
buckling will be very close to each other. Hence, the solid curve 
of k s against l /0 shown in figure 2 , faired through the data of 
table 1 , can be used to estimate k a for any value of 0 between 1 and « 


CONCUJDUKr REMARKS 


From a consideration of both symmetrical and antisymmetrical buckling 
patterns, the shear buckling stresses of rectangular flat plates with 
clamped edges have been more correctly evaluated than in previous work 
wherein only symmetrical buckling patterns were considered. Through the 
use of the LagrangLan multiplier method, which permits the computation 
of both upper and lower limits to the true buckling stress, values known 

to be within lj^ percent of the true buckling stress coefficients have 

been obtained. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., September 18, 19^7 
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APPENDIX 


COMPUTATION AIDS 


In the determinants which form the stability criteria from which 
the "buckling— etre ss coefficients are derived (reference 5)> many terms 


are 


given in the form of infinite summations; that is, / 


m=l 


where 


and 


^ 


Ajjn = 2(m£ + P 2 ^' 


This summation, for a given value of n, converges very slowly, 

approximately as however, the terms of the summation approach the 

2b r 

n 2. 

value of •? — rapidly, as m "becomes large. Thus, It may "be assumed 
that above a certain value of the Index m, say M, the terms of the 

desired series "become equal to the corresponding terms of 

Thus, for example 


” 2 

z m 

3&: 


A mn 


■ 00 M 00 

/ “^ml ~ /_ m ^ B ml + y~ 

Bt=l m=l m=Mfl 


m=l 

M 


m & 

Ami 


m=l m=l m=l 


m? 
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But 


■2 c m 2 


Than 


y j£ = y~ 

A =1 ibl 2(m 2 + P 2 ) 

= ^'(P rt + ooth Pit — Pit coth 2 pj^) 

” M M 

/ m2s ffl 1 ~ ^ + ^-(3 it + coth pit - P* coth^Jt) 

m=l m=l m=l 


In the case of P = 2, where kg ~ 10, it was only necessary to take 
M = 9 to obtain sufficient accuracy. 

All of the infinite sumnations found in the determinants for both 
the symmetrical and anti symmetrical lower— limit solutions can be 
evaluated in a similar manner. In the upper-limit determinants, there 
appear summations of a more complex type of term, such as 


0 » 

z 

m=l 



where 



h 

Am 2 


Expanding the general term gives 


m 





m 2 (m 2 + P 2 ) 2 

2(m 2 + Jm 2 ) 2 [jm 2 + 4P 2 ) 2 + lt-(m 2 + p 2 )^j 



6 


NACA TN No. 1 559 


Dividing through the numerator hy the denominator gives a sum of 
the form 


A B 

~2 If + 

mr vr 


C D 

J? " 05 + 

m m 


The first four terms of this sum are sufficient to provide a 
satisfactory approximation to the general term of the desired series 
at high values of m. Thus 



/A B C d\ 

a^M+i \ m m m m V 


M 

E 




M 

£0 


B 

T + 


m 


c d\ 



where 


and 



B B 


9.23 2 

10 • 


C 


10 


■pv _ 318.2080 6 

D a jg 


But since 


CO 



m=l 


jt 2 

T 



m=l 
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00 

7k.i 

n6 9^5 

m=l 

00 A 

y~ 

4-^J> 9*50 


tha desired summation 136002863 


2Z-4r 



2 h 

A-X- -B-X- + 

6 90 


■ -js £ _ p— a?-, 


9^5 


9450 



Similar approximations are possible for the other infinite series 
of similar form in the upper-limit determinants. 
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TABLE 1 


SHEAR- STRESS COKKETCIEKTS POR WOCCfB ASPECT RATIOS OF PLATE 


Symnetrlcal "buckling 

Anti Bywnetrical "buckling 

P 

Upper 

limit 

Lower 

limit 

Average 

MaxJjBpm 

error 

(percent) 

P 

Upper 

limit 

Lower 

limit 

Average 

Maximum 

error 

(percent) 

1.00 

*14-79 

a i4.64 

14.71 

0.54 




— 

— 

I.25 

^12.61 

°12.35 

32.48 

1.05 

— 




<«*»«*«* 

1.50 

*31-56 

G ll-.k5 

II.50 

■52 

1.50 

^.os 

e n.79 

n.93 

1.24 

2.00 

*10.60 

C 10.58 

10.59 

.10 

2.00 

f 10.36 

e 10*32 

30.34 

.20 

3-00 


C 9-57 

9.62 

■92 

3-00 

d 9-71 

e 9,64 

9.67 

.41 i 

- - -- 

00 



B 8. 9 Q 

— 

— 



e 8.98 



^ixth order determinant; p=q=3 
^Hinth order determinant} p=3^ q=Q 
^Eleventh. order determinant} p- 3 , g .«2 

%lnth order determinant; p=5, q*& 

eEleventh order determinant; p=5, q.=fi 
^Eleventh order determinant; p=7 ; 

Sottalned from exact solution In reference 




(See reference 5.) 
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